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The EtMeaP and EtMeaSb nearly triangular organic salts are distinguished from most other Pd[(dmit)2] 
based salts, as they display valence bond and no long range order, respectively. Under pressure, a supercon- 
ducting phase is revealed in EtMeaP near the boundary of valence bond order. We use slave-rotor theory with 
an enlarged unit cell to study competition between uniform and broken translational symmetry states, offering 
a theoretical framework capturing the superconducting, valence bond order, spin liquid, and metallic phases on 
an isotropic triangular lattice. Our finite temperature phase diagram manifests a remarkable resemblance to the 
phase diagram of the EtMe 3 P salt, where the re-entrant transitions of the type insulator-metal-insulator can be 
explained by an entropy difference between metal and the U(l) spin liquid. We find that the superconducting 
pairing symmetry is d ± id, and predict different temperature dependences of the specific heat between the spin 
liquid and metal. 



Quasi two-dimensional organic salts have provided a play- 
ground to study the rich physics of frustrated and strongly 
correlated systems. One fascinating example is the family 
of triangular Z[Pd(dmit)2]2 salts where Z is a monovalent 
cation and dmit refers to l,3-dithiol-2-thione-4,5-dithiolate. 
At ambient pressure, most [Pd(dmit)2]2 compounds are para- 
magnetic Mott insulators with spin susceptibility showing the 
characteristic temperature dependence of a triangular spin- 
1/2 Heisenberg antiferromagnet at high temperatures jl|]. If 
geometrical frustration is less effective, one expects long 
range order to develop at low temperatures to release the 
spin entropy. Antiferromagnetic (AF) order is a natural can- 
didate, which is indeed observed in the cases of Z= Me4P, 
Me 4 Sb, Me 4 As, Et 2 Me 2 P and Et 2 Me 2 As JH. However, no 
long range order has been detected in EtMeaSb J2J], poten- 
tially exhibiting a spin liquid phase, and in EtMeaP yl |4£| 
a spin-gapped phase with columnar valence bond (VB) order 
emerges. 

The degree of frustration plays an essential role in the va- 
riety of ground states of these materials, despite the common 
insulating behavior. In-plane spatial anisotropy, such as in the 
hopping or Heisenberg exchange integrals, is a way to con- 
trol the geometrical frustration. It turns out that the AF order 
is found in most anisotropic materials listed above, while for 
EtMesP and EtMeaSb, where AF order is absent, the lattice 
is nearly isotropic triangular til yl |5[] . 

While most [Pd(dmit)2]2 salts become metallic [4] under 
pressure, the EtMesP sa lt exhibits superconductivity with 
maximum T c appearing at the border of the VB phase seen 
in X-ray measurements [3, 4]. Re-entrant transitions of the 
form insulator-metal-insulator near the phase boundary at low 
temperatures were reported [4], where the positive curvature 
of the phase boundary indicates a large entropy contained in 
the paramagnetic insulator relative to the metallic phase. The 
complexity presented in the phase diagram of EtMeaP has 
not yet been explained within a single theoretical framework. 

In this paper, we provide a microscopic theory which 
demonstrates the emergence of the superconducting, VB or- 



dered, spin liquid, and metallic phases. We focus on an 
isotropic triangular lattice aiming to understand the phenom- 
ena observed in the EtMesP and EtMeaSb salts. In addition 
to the well-known on-site Hubbard interaction, we find a spin- 
spin interaction plays an essential role in favoring the different 
ground states reported for EtMeaP an d EtMeaSb. Our theory 
also provides a spin liquid surviving between the metal and 
VB ordered phases, superconducting pairing symmetry, and 
novel behavior of the specific heat across the metal-insulator 
transition. 

A common theoretical approach for Mott insulators is to 
start from a localized Heisenberg spin model, which is ap- 
propriate for insulators. However, EtMeaP becomes metal- 
lic/superconducting when pressure is applied, exhibiting a 
metal-insulator (MI) transition. Capturing the MI transition 
requires an electronic description, such as a Hubbard model 
at an intermediate coupling |g-|§|]- Can one access the com- 
petition between spin liquid and VB ordered insulators or be- 
tween metal and superconductor within the Hubbard model 
near the MI transition? Prior work IICJ 1 1111 suggests diffi- 
culty in achieving a superconducting phase within the pure 
Hubbard model at intermediate coupling in contrast to when a 
spin-spin interaction is included. Combining the experimen- 
tal facts and the previous theoretical work, we suggest the fol- 
lowing t-J-U Hamiltonian is a minimal model for isotropic 
triangular [Pd(dmit)2]2 salts. 
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where J2iii) runs over nearest neighbors, t > is the near- 
est neighbor hopping (we will set t = 1 as our energy scale 
below), U is the on-site Hubbard repulsion, and J is the spin- 
spin interaction. We emphasize that the J term is being con- 
sidered as an independent effective interaction. The ci a and 
c ia . are the electron operators and the electron spin operator is 
Si = ^2 a » cl a (&)apCii3, where a are the Pauli matrices. 

We take a slave rotor approach first developed by Florens 
and Georges J12lll3n which can handle the intermediate cou- 
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FIG. 1 : (a) The doubled unit cell (b) Form factors used in the mean 
field equations. 



pling regime where charge fluctuations are important. In the 
slave-rotor representation [13], the electron operators are ex- 
pressed as Ci a — fi a e t6i , where exp (—i9.j) is the rotor lower- 
ing operator and fi a , f* ia are fermionic spinon operators. Us- 
ing this representation the Hamiltonian (fTJ is written as 111 311 



H = -t £ flf ja e-^-^ 
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where Sj = 1/2 Y, a p fia(°)<*pfiP is the spinon spin opera- 
tor and Li is the rotor charge operator 11311 . It is important that 
the constraint of Li + ^ Q flafia = 1 be satisfied to project 
the Hilbert space into the physical subspace at half-filling. 

We adopt the slave rotor theory and extend to a enlarged, 
two site unit cell to study the competition between uniform 
states and states which minimally break physical translational 
invariance. The spinon operators for the two sites in our 
unit cell are denoted by a a and di^ a as indicated in Fig. 
|l(a)| Introducing the mean fields, Xij = |E,(/l/j»). 
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particle-hole Nambu basis. The additive constant is the Fermi 
sea energy. Treating the rotors as a bosonic variable, i.e., 



oi&i 



bi subject to the constraint of 



1, the rotor 



part of the Hamiltonian can be represented by the (imaginary 



time) action, 
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where the final form is obtained as in the spinon case by de- 
composing into the two sublattices and Fourier transforming. 
The iXi(= h) enforces the constraint and acts as a chemical 
potential. To consider condensation, we separate out the low- 
est energy component of the lower band 7^7 by diagonalizing 
the above action, and assuming k = condensation can be 
treated as an independent variable. Denoting the condensate 
density as \z\ 2 , the saddle point conditions for h a , hi, and z 
give constraint equations for the number of bosons. 
The full set of mean field equations is given as 
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where a = ±, the ff"(fc) form factors with r = 1 to 
6 are defined in Fig. |l(b)| E[(k) are the eigenvalues of 
Mk and Ut is the matrix of eigenvectors. The 4x4 ma- 
trices Ox(k) are defined by the schematic relation X = 
~k Efe(^[^xw^fe). The bosonic distribution rn,{x) is de- 
fined as n;, (a;) = J t£- coth \fl\ ^r ) and nf (x) is the Fermi 
function. The constraint equations for the number of fermions 
and bosons are given by 
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where s = c,d, the form factors uf(k) are the components 
of the rows of the transformation matrix that takes the rotor 
action to diagonal form and n c , rid are the densities on each 
sublattice. When the condensation amplitude is finite we also 
have Jq\z\ = 0. 

Fig. [2]shows zero temperature phase diagram varying J and 
U obtained by solving the above coupled mean field equations 
Eq. |2] self-consistently while enforcing the constraints, Eq. 
[3] This involves 12 complex mean fields, 6 \ij an d 6 By as 
indicated in Fig. |l(a)| linking the fermionic and bosonic mean 
field theories of the spinon and rotor sectors respectively. For 
the pairing we use an ansatz of A r = | A\e l6s ^ with uniform 
magnitude but allow the phase to vary within the unit cell. 
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FIG. 2: (Color online) The J vs U phase diagram. The dotted line 
denotes a second-order transition, while the solid line is first-order. 
The label d + id denotes the superconducting phase, and the inset 
shows the phase of A on each bond. The (7(1) spin liquid is labelled 
SL. The pattern of \ an d B in the dimer state is inset in the dimer 
phase, with thick lines denoting large values and thin lines small val- 
ues. We have rescaled U — > U/2 as discussed in 11311 . 



Four phases are identified in Fig. [2] a superconductor, a 
metal, a spin liquid and the dimer state. The uniform metal- 
lic state is characterized by the condensation of the bosonic 
degree of freedom and a uniform, real Xij — Ixl implying 
the existence of an electron Fermi surface and gapless charge 
excitations. On the insulating side, the spin liquid is also char- 
acterized by a uniform, real Xij but g a PP e d bosons, meaning 
the existence of a spinon Fermi surface but gapped charge ex- 
citations. The dimerized phase is characterized by \i = 1/2 
and other \i = (with a similar structure for B) and A = 0. 
Since the mean field solution consists of completely isolated 
bonds, any dimer covering of the triangular lattice will have 
the same mean field energy at this level of approximation. In 
the superconductor we find the pairing is d ± id, as the phase 
winds by 27r/3 around a triangle. It is fully gapped, with no 
nodes. Note that we do not find a Z2 spin liquid, that is a spin 
liquid phase with finite spinon pairing but uncondensed rotor 
bosons. 

One can show via perturbation theory that if U/t 3> 1, 
Bij ~ Atx*j/U, so the spinon Hamiltonian becomes iden- 
tical to a slave-fermion description of the Heisenberg model 
with J = 4t 2 /U. The result of Rokhsar [ 14] states that under 
our conditions a dimerized configuration is among the mean- 
field ground states. Thus we expect a dimerized state at large 
U. However, the U(l) spin liquid survives between the metal 
and dimer phases near MI transition due to charge fluctua- 
tions. The window of the spin liquid phase is quickly wiped 
out as we increase J. One expects that as J becomes large 
our mean field solution will be equivalent to that of a Heisen- 
berg model decoupled in \ ar, d A channels. The argument of 
Rokhsar [14] can be generalized to this case, and implies that 
the mean field ground states include dimerized states, where 
Xij = A,j = on some bonds while on others \xij\ 2 + \^ij | 2 



is maximal. Thus we expect that the dimerized state will be 
among the ground states as we increase J. Lee and Lee [15] 
have looked at the pure Hubbard model ( J=0 in our model) 
within the same mean field framework discussed above. How- 
ever, they considered only translationally invariant solutions, 
and found a continuous transition from a metal into a spin liq- 
uid followed by a first order transition into a 7r-flux state. We 
find the 7r-flux state is replaced by the dimer phase, and the 
window of spin liquid is also narrower than that of Ref.[ 15]. 

It is straightforward to compute finite temperature phase 
diagram as shown in Fig. [3] (a) and (b). One peculiarity is 
that there is no condensation of the bosons in two dimensions 
at any finite temperature, and thus no metallic phase. This 
is an artifact of the spherical approximation, that is treating 
the rotors as bosons 1 1311 . As the rotor action is equivalent to 
one used to describe Josephson junction arrays, we expect a 
Kosterlitz-Thouless type transition from metal to spin liquid, 
with quasi-long range order at low temperatures. In fact, in- 
cluding gauge fluctuations, the nature of the phase transition 
from metal to spin liquid may be a crossover Ill6lll7ll . In Fig. 
(a) and (b) we include a weak coupling between the lay- 
ers via — 2t z cosk z in the dispersion with t z ~ 10~ 4 , which 
should be present in quasi two-dimensional organic salts, to 
provide a consistent bosonic description for the whole phase 
space. 

The J = 0.25 phase diagram in Fig. [3] (a) show a strong 
resemblance to the phase diagram of EtMcaP found in Refs. 
[4, 1 8] . Identifying the VB ordered state with the dimer state 
deserves some discussion. First, since the dimer states are 
completely degenerate, any small anisotropy should break the 
degeneracy and pick out a single dimer configuration, result- 
ing in the VB ordered phase seen in Ref. J2B. A coupling 
to the lattice should play an important role in this context as 
discussed in [4]. The spin susceptibility shows the expected 
behavior as well, being exponentially suppressed in the dimer 
phase. Second, the T^ 1 data for the VB ordered state indi- 
cates a significant amount of inhomogeneity, possibly due to 
disorder. Taking this disorder into account could drive the 
first order transition we have found theoretically to the sec- 
ond order transition reported in [4]. As well, one expects that 
including corrections beyond our mean field theory could be 
qualitatively similar to a quantum dimer model Ill9ll20ll . It has 
been argued [21] that a VB solid phase in these models could 
show a second order transition to a VB liquid phase at finite 
temperature, analogous to our dimer to spin liquid transition 
becoming second order. At higher temperatures, the spin liq- 
uid reappears because of the large entropy associated with the 
spinon Fermi surface and the gapped nature of the dimer state. 
Note that near U ~ 4.2—4.6 as temperature increases the tran- 
sitions, insulator (dimer) to metal to insulator (spin liquid) oc- 
curs, similar to the re-entrance seen by resistivity experiments 

At lower values of J the spin liquid state persists even at 
zero-temperature, a scenario which has been discussed ex- 
tensively for the organic compound k — (ET)2Cii2(CN)3 
IU5J l22l 12311 . The J = finite temperature phase diagram 
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FIG. 3: (Color online) Finite temperature phase diagrams, with 
J = 0.25 for (a) and J = 0.00 for (b). Curves have been fitted to 
the boundaries to guide the eye. A dotted line denotes a second-order 
transition, a solid line is first-order and a dashed line is a crossover. 
Note the superconducting phase at small U for J = 0.25 in (a). In 
(c) specific heat as a function of T is shown for various values of U 
denoted by different colors in (b). Note that the slope of the curves 
increases with U indicating more density of states in spin liquid in- 
sulator than metal. 



120° degree magnetic ordering. Although they treated the 
rotor sector differently, at large U/t both methods give the 
same asymptotic relation By ~ At%ij* /U, which allows us 
to estimate the critical U/t « 18, above which the magnetic 
order sets in for J/t = 0.25. Since we are only interested 
in the intermediate U, our main results near the MI transi- 
tion are not affected by magnetic ordering of this type. Note 
that the tendency to order magnetically is also affected by the 
anisotropy as we discussed above, which distinguishes EtMe.^ 
P and EtMe3Sb from other [Pd(dmit)2]2 based organic salts. 

In summary, we propose an effective model for the isotropic 
triangular Pd(dmit)2 salts. We show that the Hubbard U and 
spin exchange interactions account for the diversity present in 
the phase diagram of the EtMeaP salt, specifically the prox- 
imity of superconductivity and VB order near the pressure in- 
duced MI transition. The spin liquid is realized for weak spin- 
spin interaction, which may be relevant for the EtMeaSb salt. 
Our finite temperature phase diagram shows the re-appearance 
of the spin liquid at finite temperatures due to entropic rea- 
sons, explaining the re-entrant insulator-metal-insulator tran- 
sitions. We also present the temperature dependence of the 
specific heat, which is qualitatively different between the spin 
liquid and metal and can be tested in future experiments. 

We thank M. Hermele, Y.B. Kim, A. Paramekanti and S.S. 
Lee for useful discussions. This research was supported by 
NSERC of Canada and the Canada Research Chair program. 



is shown in Fig. [3(b) and has a window of spin liquid be- 
tween a metallic and dimer phase, which gets wider at finite 
temperature due to the larger entropy contained in the spin liq- 
uid than the dimer phase. It is tempting to argue that EtMeaSb 
falls into the spin liquid regime, while EtMeaP belongs to the 
dimer phase nearby the spin liquid due to different exchange 
strengths J. As we discussed, the gauge fluctuations are ex- 
pected to turn the metal-spin liquid transition into a crossover 
rather than a sharp transition, but the details are still under de- 
bate. Due to the paramagnetic state appearing above the VB 
ordered phase at temperatures comparable to this crossover 
scale, reliable information on the gauge theory is needed to 
extract any experimentally relevant information. However, we 
find that even at the level of mean field theory, there is differ- 
ence in temperature dependence of the density of states in spin 
liquid and metal phases. The specific heat is plotted in spin 
liquid and metal in Fig. 0(c). We note the linear behavior in 
both metal and insulating phase, with increasing slope as we 
pass from metal into the spin liquid. This indicates a larger 
density of states at the Fermi level in the insulator, in con- 
trast to conventional MI transition where the density of states 
decreases to zero. 

It is worthwhile to comment the absence of magnetic order- 
ing in the phase space considered here. At large U or J one 
expects magnetic ordering of 120° rotated spins. The mag- 
netic order was considered in [24] where the spin-spin term 
is decoupled into the \ an d magnetic channels, imposing a 
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